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' In this paper, we propose an efficient scheme to drive three atoms in an optical 

cavity into singlet state via adiabatic passage. Appropriate Rabi frequencies of the 
classical fields are selected to realize present scheme. The scheme is robust against 
the deviations in the pulse delay and laser intensity through some simple analysis of 
adiabatic condition. It is notable that the estimated range of effective adiabaticity 



a* 



qh. condition coincides with the numerical results. When taking dissipation into account, 

we show that the process is immune to atomic spontaneous emission as the atomic 
excited states are never populated in adiabatic evolution. Moreover, under certain 
^ ■ conditions, the cavity decay also can be efficiently suppressed. 

o 
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I. INTRODUCTION 

The technique of adiabatic passage [l-3] has been proved as an effective coherently con- 
trol dynamical process to realize quantum information processing (QIP) in many schemes 
4-llj]. Unlike the stimulated Raman adiabatic passage (STIRAP) where the Stokes pulse 
vanishes first, the two pulses vanish simultaneously while maintaining a constant finite ratio 



of amplitudes in f-STIRAP |12Hl4|. which guarantees the creation of any pre-selected co- 
herent superposition of ground states. A remarkable superiority of STIRAP and f-STIRAP 
is that if the evolution is adiabatic, the states of the system evolve within an adiabatic dark 
state subspace and the relevant states contain no contribution of the excited atomic states, 
so the spontaneous emission from excited states can be suppressed. Another advantage of 
adiabatic passage is the simpleness, for it needs not consider the precise tuning of pulse 
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areas, pulse widths, pulse shapes, pulse delay and detunings. 

Cabello has proposed a special type of entangled state which is called iV-particle iV-level 
singlet states in 2002 15]. These states are key solutions to many problems such as N- 
strangers, secret sharing and liar detection, which have no classical solutions. It has been 
shown that some types of supersinglets may violate the local hidden theory 16] and can 
be used to construct decoherence-free subspaces (DFSs), which are robust against collective 
decoherence |17|. Nevertheless, the generation of these states for N = 3 which have the form 

4=(|012) - |102) - |210) + |120) + |201) - |021», (1) 

is still a problem both in theory and experiment. Here |0), |1) and |2) represent three ground 
states. Jin et al. have proposed a scheme of generating a supersinglet of three three-level 
atoms in microwave cavity QED based on the resonant atom-cavity interaction isj]. The 
scheme is sensitive to the loss of cavity since the three atoms are sequentially sent through 
three different cavities and cavity fields which act as memories. Lin et al. also have raised a 

n 

protocol for the preparation of a singlet state with three atoms via Raman transitions [19]. 
However, there will be a considerable influence caused by the cavity decay and spontaneous 
emission of the atoms. Shao et al. have put forward an approach by converting two-atom 
singlet state into three-atom singlet state via quantum Zeno dynamics 20], which needs to 
control the interaction time exactly. 



To overcome these problems in Refs. 



1 81420]. we propose a scheme where the state of the 



system evolves within a dark-state subspace via adiabatic passage. Three Gaussian shape 
pulses are used to complete the scheme and by this way the two Rabi frequencies of the laser 
fields maintain a constant finite ratio of amplitudes to get the desired final singlet states. 
Also sufficient adiabaticity can be achieved by choosing such pulses, for the estimates based 
on simple analysis of adiabaticity condition show the scheme is robust against deviations 
such as the pulse delay and pulse intensity, and the analysis is also validated by numerical 
calculation. Moreover, we also analyze the influence of the choice of laser intensity on the 
dissipation and achieve a relative high fidelity by choosing appropriate parameters. 

The paper is organized as follows. In section II, we introduce the models for adiabatic 
passage and the pulses to realize the scheme. In section III, we discuss the robustness of 
parameters mismatch and dissipation due to the atomic spontaneous emission and cavity 
delay. Section IV contains the concluding remarks. 
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II. GENERATION OF THREE-ATOM SINGLET STATE 

As depicted in Fig. 1, three four-level atoms with tripod configuration are trapped in a 
bimodal vacuum cavity field. Each atom has an excited state |e) and three ground states 
|/r) and \r). Supposed that the transition between the levels \e)i -H- |/z,)i(|/ii)i) (i = 
1, 2, 3) is resonantly coupled to the cavity mode with the coupling strength g^gui) and the 
transition |e)j -H- |r)j is resonantly driven by the classical pulse with the Rabi frequency 
In the interaction picture, the Hamiltonian for the whole system can be written as {h = 1) 

H tot = Hi + H c , 



Hi = ^2ni(\e)i(r\ + \r)i{e\), 

i=l 
3 

H c = ^giL(a L \e)i(f L \+ a[\f L )i(e\) + g iR (a R \e}i(f R \+ a^ R \f R )i{e\), (2) 



i=3 



where a\{a R ) and dL{a R ) are the creation and annihilation operations for the left(right)- 
circular polarization cavity mode, respectively. We assuming g^ = g iR = g, fi 2 (t) = ^(t) = 
fl(t), and fli(t) to be real in the present paper for simplicity, the whole system will evolve 
in the following closed subspace: 



0i) = |r/ L /fl)i 23 |00) \<h) = |r/B/L)l23|00) 



i 3 > = |/tr/B>i 23 |00> 



a L a R i 



5 ) = |/L/Rr) 123 |00) 



h) = |e/ L /fl)i23|00) 



a h a R i 



09) = \f L ef R )l23\00) aL a R , |0io) = | iW L ) 123 1 00) 



/a L a R i 



k) = \f R rf L )i23\00) 



a L a R i 



f>e) = \f R hr)i23\00) 



O.LO.R 1 



8 ) = |e/fl/L>i23|00) 



a L a R i 



la L a R i 



1012) = |/B/Le)i23|00) 



a-LO-Ri 



013) = \fLfLfR)l23\±0)a L a R , |0u) = | /l/r/l) 123 1 10) a L a R , 
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'is) = |/b/z,/z,)i28|10)o £ ok, I0ie) = |/b/r/l)i23|01) 



a La r , 



1017) = \fRfLfB)l23\0l)a L a R , |0is) = \f if rJr) I2l\®l) a L a R , (3) 

the subscripts 1, 2, 3, cll and an represent atom 1, atom 2, atom 3, left-circular cavity mode 
and right-circular cavity mode, respectively. 

There are six dark states with null eigenvalue in this subspace. We orthogonalize these 
states and get a special dark state \S) which will evolve into an independent subspace while 
other states remain unchanged. This state can be expressed as 

\S}= ^ [ ~ mg(\h) - 102}) - ^|^(|03) - |04) - |0 5 ) + 106)) 

+^^(|0i3>- |0u)- |0i 6 ) + |0ir))l, (4) 



where N = yj£l\{t)g 2 + 2Q?(t)g 2 + VL\(t)VL 2 (t). Note that the state in Eq. (4) contains no 



contribution from the atomic excited states 



21 



22] , which can be considered as unpopulated 



during the whole interaction process, if the evolution is adiabatic. When the system is 
initially in the state 

10*} = ^(I0i) - 102}), (5) 

under the condition 

^>fi!(t),n(t), (6) 

and the Rabi frequencies following such behaviour 

hm — , = 0, hm — . = tana, 7 

we have an approximate evolution process from the initially entangled state |0j) to three- 
atom singlet state 

|0t) = ^(I0l) - 102) + 103) - 104) - |0 5 ) + 106)), (8) 

when a = arctan2, which is the result we need. Let us define that 9{t) = arctan[f2i(t)/fi(£)], 
note that the rate of the change of the mixing angle 9{t) must be much smaller compared to 



the smallest separation Au(t) of the corresponding eigenvalues the specific expression 

is 

\6{t)\ < Aco{t). (9) 

Under these conditions, the evolution is adiabatic, and the system will remain in the dark 
state \S). 

Next we discuss the Rabi frequencies and other correlative parameters that make the 
scheme experimentally feasible. We use three time-dependent pulses as follows 

Qx{t) = sinaft exp-(*- r ) 2 / T2 , 

Q(t) = cosafioexp"^ 2 ^ 2 + ft exp"^ 2 ^ 2 , (10) 

and the variation of the two time-dependent Rabi frequencies fl\ (red dashed curve) and Q 
(blue solid curve) of lasers for atoms is shown in Fig. 2(a). By choosing r = QO/g, T = 
80/ g, Qq = 0.2 g, we can see that when t > 100/ g, the two Rabi frequencies are approximate 
to meet the relation Qi ~ ^Q. Thus the initial state \4>i) will transfer into the target state 
\4>tj- Fig. 2(b) shows the time evolution of the populations of the components of the state 
\S) P\(P2) (green dashed curve), P 3 (P4, P 5 , P 6 ) (purple dotted curve) and Pi 3 (P 14: , P w , Pi 7 ) 
(red solid curve), respectively. 

It is obvious that the population is almost completely transferred from the state \<pi) to 
the state \(f> t ) without populating other states during the dynamical process, which means 
the influence of the atomic spontaneous emission is effectively suppressed. For the density 



matrix of the two entangle states pi, p2, the Bures fidelity can be defined as [23] 
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F=(tr^p 1 1 /2 p 2 p\ /2 ) , (11) 

The relation between the fidelity and the evolution time t is shown in Fig. 2(c). When 
T = 100/g, the fidelity of the three-single state is 0.9996, which means we finally realize 
an almost perfect target state. Note that we need not control the time accurately since a 
fidelity F higher than 0.8775 can be obtain when t > 20/ g. 

III. ADIABATICITY CONDITION AND NUMERICAL ANALYSIS 



Now we come to discuss the conditions that a system with the initial state \<f>j) evolves 
adiabatically into the target state \4>t)- in our analysis of the robustness of adiabatic passage 
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against variations in the experimental parameters, we start with the adiabatic condition (9). 
For the restrictive conditions of Eq. (6) and the given pulse shapes we have 

4r ^(t)sino; 



0(t) = 



T 2 sin 2 a + [cos a + £(t)] 2 ' 



' " w 2 2 2 2 



fi 2 (t) 2fi 2 (t) 



= n e'^y T2 f^ + l[cosa + m] 2 = ^f f (t), (12) 

where = e - 4rt / T2 . Then a limit for the pulse delay can be get by analysing the upper 
expressions. 

We know that non- adiabatic conditions are most likely to occur when 9(t) reaches a 
maximum, that is £(£o) — 1- Take t = for example, then Eq. (12) is equal to 

2r 1 

0(*o) = 0m*r(*) = ^ tan-a, 



^ / x ^ ^2 / sin 2 ct 2 . . _ . . 

n eff (to) = n oe - T y ^— + 3 (cos a + 1) 2 . (13) 

The maximum of 9 (t) increases with r (see Eq. (13)), while fi e //(^o) is not necessarily the 
maximum of Q e ff(t). To reach the adiabatic conditions (9) we must have Q e ff(to) > n9(to), 
where n is a 'sufficiently large' number and the choice of which depends on how much 
non-adiabaticity can be allowed. Then we find an upper bound on r, 

n r > 2nr - tan ^ ^ (14) 

T y^ + Kcosa + 1) 2 

We can see that the Rabi frequency needs to increase exponentially with r to suppress the 
non-adiabatic transitions. 

Note that both 9(t) and Q(t) are pulse-shaped, it is convenient to find their full widths 
at half maximum(FWHM) for the following discussions. 



T 2 / / 1 1 
Tg ~ — In I W 1 + cos 2 -a + cos -a 
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T n ~ 2r + 2Tv / hi2. (15) 

When r — >■ 0, the pulses of will be broaden and get broader than Q(t) . Then condition 
(9) will be violated in the early time and late time during the evolution. This problem can 
be solved by controlling the width of 9 (t) smaller than the width of Q(t), that is Tg < T n . 
Then we obtain a lower bound for r. When a = arctan2 it reads as 

r > 0.25T. (16) 

From the above analyses we can see that although for any sufficiently strong laser pulse 
and delay r > the f-STIRAP should work, there exists a more accurate adiabatic condition. 
For n = 5, the range is 0.25T < r < 0.97T for Q T = 12, 0.25T < r < 1.06T for Q T = 16 
and 0.25T < r < 1.13T for tt T = 20. We also plot the relationship of the fidelity F versus 
the ratio r/T in the situation that VL$T = 12 (blue dotted curve), QqT = 16 (green solid 
curve) and QqT = 20 (red dashed curve) by solving the master equation numerically in Fig. 
3(a). The figure coincides with the results we deduce from Eq. (9) almost perfectly, which 
verifies that the adiabaticity is most easily achieved under the range we obtain. Fig. 3(a) 
also implies that for a particular value of fi T, there exists a relative wide range of r/T, 
which means that our protocol is robust against the pulse delay in practical experiment. 
Fig. 3(b) shows the change of the time dependence of the smallest separation Au(t) (blue 
dotted curve), the simplified smallest separation Q e ff(t) (green solid curve) and the mixing 
angle 9{t) (red dashed curve) with the chosen parameters, which implies that the parameters 
our scheme chooses fit condition (9) well. 

In all the above discussions, we have not considered any dissipation and assume the 
system does not interact with the environment. However, the system will interact with the 
environment inevitably which has influence on the availability of our scheme. Hence we will 
focus on discussing the influence of dissipation induced by the atomic spontaneous emission 
and the cavity decay. When we consider decoherence, the master equation of motion for the 
density matrix of the whole system can be expressed as 

p = -i[H tot , p] - ^{a\a L p - 2a L pa[ + pa[a L ) 
~ -rr( a R a RP ~ 1a R pa ] R + pa ] R a R ) 
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3 ^ -pfc 

k=l m=f L ,g,f R 

where T^ m is the spontaneous emission rate from the excited state |e) to the ground states 
\m) (m = /z,,r, /r) of the fcth atom, kl(kr) is the decay rate of the left(right)-circular 
cavity mode. We assume T^ m = T = T /3 and Kl = Kr = k for simplicity. Fig. 4 shows the 
relationships of the fidelity F versus the ratios Qo/g and T/g, VL Q /g and n/g, respectively. 
The spontaneous emission rate makes a slighter influence on F under a larger laser intensity 
while cavity decay causes an opposite situation. The physical mechanism behind the behave 
of the former case is just the adiabatic condition (9), which is satisfied better along with laser 
intensity increasing. That means under a relative large laser intensity, i.e. at flo/g — 0.3 in 
Fig. 4., the passage is more likely to evolve within the adiabatic dark state subspace which 
does not involve the excited atomic state |e). Therefore, the dissipation caused by atomic 
spontaneous emission decreases. However the populations of states where the cavity field 
are excited increase with laser intensity according to Eq. (4), which increase the dissipation 
caused by cavity decay finally. In such a way, the change of laser intensity decreases one error 
source while increasing another. So an appropriate value Qq should be chosen when taking 
both the two factors into account. We also plot the relationship of the fidelity F versus the 
ratios n/g and T/g by solving the master equation numerically in Fig. 5. Therefore we can 
see that under certain conditions both atomic spontaneous emission and cavity decay have 
a slight influence in fidelity F, since for a large atomic spontaneous emission T/g = 0.05 and 
cavity decay n/g = 0.05, the fidelity is still about 0.9244. Therefore our scheme is robust 
against the two error sources and achieve a superior result in theory. 

Finally, we give a brief discussion about the basic factors for the experimental realiza- 
tion. The atomic configuration might be achieved in cesium atoms in our scheme. The 
state \r) corresponds to F = 4, m = 3 hyperfine state of 6 2 Si/ 2 electronic ground state, \fi) 
corresponds to F = 3,m = 2 hyperfine state of 6 2 Si/ 2 electronic ground state, |/#) corre- 
spond to F = 3,m = 4 hyperfine state of 6 2 5*i/2 electronic ground state, |e) corresponds to 
F = 4, m = 3 hyperfine state of Q 2 P\/2 electronic state, respectively. In recent experimental 



condition 



24 



251, the parameters g = 2ttx750M Hz, T = 2nx2.62M Hz, re = 2nx3.5MHz 



and the optical cavity mode wavelength in the range between 630 ~ 850 nm is predicted to 
achieve. By substituting the ratios n/g = 0.0047, T/g = 0.0035 into Eq. (17), we will obtain 
a high fidelity of about 0.956, which shows our scheme to prepare three-atom singlet state 
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\4>t) is relatively robust against a realized one. 

IV. CONCLUSION 

In summary, we have proposed a scheme to generate a three-single state for three atoms 
trapped in an optical cavity via the adiabatic passage of dark state. The significant feature 
is that we need not to control the laser time exactly and it is robust against variations in the 
laser parameters such as pulse delay and laser intensity. So the scheme is robust, effective 
and simple. When considering dissipation, we can find that the protocol is robust against 
atomic spontaneous emission since the states evolve in a closed subspace where the atoms 
remain in the ground states in a adiabatic evolution. By choosing proper parameters, the 
scheme is also insensitive to cavity decay by numerical calculation intuitionally. The result 
shows that the scheme have a high fidelity and may be possible to implemented with the 
current experiment technology. 
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FIG. 1: The level configuration of the scheme. The transition |e) — > \fi) and |e) — > |/#) are 
coupled to left-circularly and right-circularly polarised cavity modes, respectively. A classical laser 
driver the transitions |r) — > |e). 
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FIG. 2: (a) The time dependence of the laser fields for atoms. Here the red dashed curve and the 
blue solid curve represent the Rabi frequencies Oi and f2, respectively, (b) Time evolution of the 
populations. Here the green dashed curve, the purple dotted curve and the red solid curve represent 
the populations P\{P2)-, Pz[P&-,P$-,Pg) and Pi3(Pu, Pw, Pn), respectively, (c) Time evolution of 
the fidelity. We have chosen r = 60/g, T = 80/g, W 2 /g = 100, ft = 0.2g. 
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FIG. 3: (a) The fidelity F vs. the ratio t/T when OoT is in different values. Here the blue 
dotted curve, the green solid curve and the red dashed curve represent the values f2oT = 12, 16 
and 20, respectively, (b) The time dependence of the the change of the smallest separation Aui, 
the simplified smallest separation f2 e // an d the mixing angle 9. Here the blue dotted curve, the 
green solid curve and the red dashed curve represent the change of Aw, ^ e // an d 9, respectively. 
We have chosen r/g = 60, W/g = 80, n /g = 0.2. 
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FIG. 5: The infuence of ratios n/g and T/g on the Fidelity F of the three-atom singlet state. 



